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Abstract 

In this paper, we present two classes of lopsided systems and discuss their ana¬ 
lytic integrability. The analytic integrable conditions are obtained by using the 
method of inverse integrating factor and theory of rotated vector field. For the 
first class of systems, we show that there are n + 4 small-amplitude limit cycles 
enclosing the origin of the systems for n > 2, and 10 limit cycles for n = 1. 
For the second class of systems, we prove that there exist n + 4 small-amplitude 
limit cycles around the origin of the systems for n > 2, and 9 limit cycles for 
71 = 1. 

Keywords: Nilpotent Poincare systems; analytic integrability; Lyapunov 
constant; Rotated vector field. 


1. Introduction 

Integrability is one of the most important and difficult problems in studying 
ordinary differential systems. To explain the problem, consider a planar analytic 
differential system, described by 

” = yi AA’ (in 

v = u + V (u, v), 

where dot indicates differentiation with respect to time t , U and V are real 
analytic functions whose series expansions in a neighborhood of the origin start 
at least from second-order terms. By the Poincare-Lyapunov theorem, system 
HU has a center at the origin if and only if there exists a first integral, given 
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in the form of 


( 1 . 2 ) 


4>(u,v) = U 2 + V 2 + ^2 4>kjU k V j , 
fc+j=3 

where the series converges in a neighborhood of the origin. Determining whether 
the origin of system ED is a center or focus is called center problem. Another 
important problem in study of system ED is the existence of analytical first 
integral in a small neighborhood of the origin of system ED- If there exists such 
an analytical first integral, the origin of system ED is a center, in particular, 
called an analytic center, see [lj. 

It is well known that it is difficult to distinguish focus from center when the 
singular point is degenerate. Many research works have been done in this di¬ 
rection. For example, analytic systems having a nilpotent singular point at the 
origin were studied by Andreev [ 2 ( in order to obtain their local phase portraits. 
However, Andreev’s results do not distinguish focus from center. Takens Q] 
provided a normal form for nilpotent center of foci. Later, Moussu jd| found the 
C°° normal form for analytic nilpotent centers. Further, Berthier and Moussu 
[fj studied the reversibility of nilpotent centers. Teixeria and Yang [fi] analysed 
the relationship between reversibility and the center-focus problem, expressed 
in a convenient normal form, and studied the reversibility of certain types of 
polynomial vector fields. Han et al. considered polynomial Hamiltonian sys¬ 
tems with a nilpotent singular point, and they obtained necessary and sufficient 
conditions for quadratic and cubic Hamiltonian systems with a nilpotent sin¬ 
gular point which may be a center, a cusp or a saddle, see In particular, 
the local analytic integrability for nilpotent centers was investigated [ 8 ( , for the 
differential systems in the form of 

x = y + P 3 (x,y), 

y = Q3(x,y) 7 

which has a local analytic first integral, where P 3 and Q 3 represent homoge¬ 
neous polynomials of degree three. For third-order nilpotent singular points of 
a planar dynamical system, the analytic center problem was solved by using the 
integrating factor method, see for example 

The Kukles system, as a well-known example, has been investigated inten¬ 
sively on the existence of its limit cycles as well as its integrability. For the 
following particular Kukles system, 


x = y, 

y = —x + a 3 x 2 + a 2 xy + a 3 y 2 + 04 a ; 3 + a 3 x 2 y + a§xy 2 + CI7?/ 3 , 


the conditions under which the origin of the system is a center have been exam¬ 
ined in EE EH HE El. More details about the Kukles system can 
be found in jl7l |. The so-called extended Kukles system, 


x = y(l + kx), 

y = —x + a\X 2 + a-ixy + a 3 y 2 + a 4 X 3 + a 3 x 2 y + a^xy 2 + a-?y 3 , 
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has also been considered to obtain the center conditions Bill- Recently, 
center problem for some more generalized Kukles type systems have been studied 
[2(1 [2lt I22I ]. A kind of Lienard systems of type (n, 4) for 3 < n < 27 was 
investigated and they obtained the lower bound of the maximal number of limit 
cycles for this kind of system in [231. 

Research on Hilberts sixteenth problem in general usually proceeds by the 
investigation on specific classes of polynomial systems, much effort has been 
devoted in recent years to the investigation of various systems such as poincare 
system, Able equation, lopsided system and so on. The Kukles system is perhaps 
the earliest example of lopsided systems which have the following forms 


or 


x = -y, 

y = x + P{x,y), 

x = -y + P(x,y ), 

V = x. 

Since then, lopsided systems have drawn more and more attention to researchers. 
Lopsided quartic and quintic polynomial vector fields have been studied and 
center conditions were obtained [24], 251. Furthermore, Gine (Hj proved that 


there is exactly one isochronous system for lopsided quartic system, and the 
origin never can be an isochronous center for lopsided quintic system. For 
seventh-degree lopsided system Soriano and Salih [27j showed that the origin is 
a center if and only if the system is tinre-reversible and if it is not, no more than 
seven local limit cycles can bifurcate from the origin under certain conditions. 
However when the origin is a degenerate singular point, there are fewer results 
because it is difficult to compute the Lyapunov constants. The cubic lopsided 
system with a nilpotent singular point has been investigated intensively. For 
example, Alvarez and Gasull [28| proved that three limit cycles can bifurcate 
from a nilpotent singular point of the following system: 


y = a\x 2 + a, 2 xy + a%y 2 + 04 a ; 3 + a$x 2 y + a^xy 2 + ary 3 , 

via an analysis based on normal forms. Then, Liu and Li [29j showed that by 
making a small perturbation to the linear terms of it can exhibit four 

small-amplitude limit cycles. Bifurcation of limit cycles and center conditions 
for the following two families of lopsided systems with nilpotent singularities, 

x = -y + P 4 {x,y), 
y = -2a; 3 , 

and 

x = -y + P 5 (x,y), 
y = -2a; 3 , 

have been considered by Li et al. (30!> where P±{x,y) and P$(x,y) represent 
homogeneous polynomials in x and y of degree four and five, respectively. Their 


3 






results show that it is more difficult to distinguish focus from center when the 
singular point is degenerate. As far as analytic center of lopsided system is 
concerned, it is more challenging to distinguish it from focus. So, in this paper, 
we shall discuss analytic center conditions and bifurcation of limit cycles for two 
classes of lopsided systems with a cubic-order nilpotent singular point, given by 

x = y + H 3 (x,y) + H 2n+ 3 (x,y), 

V = -2x 3 , 1 ' j 


and 

± = V ' (1.5) 

y = —2x 3 + H 3 (x, y) + H 2n + 3 ( 2 , V ), 

where Hk(x,y) represent a /cth-degree homogeneous polynomial in x and y. 

The main goal of this paper is to apply the method of integrating factor and 
theory of rotating vector fields to distinguish analytic integrability conditions 
and to find the conditions for analytic centers. This work is a continuation of 
that for the Kukles system with a degenerate singular point. In next section, 
we present some known results which are necessary for proving the main result. 
We derive the analytic center conditions for the centers of systems m and 
m in Sections 3 and 4, respectively. Finally, conclusion is drawn in Section 
5. 


2. Preliminary results 


In this section, we present some relative notions and results taken from |3ll |32|, 
which will be used in the following sections. A system whose origin is a cubic- 
order monodromic singular point can be written as 


x = y + px 2 + <HjX l y 3 = X(x,y), 

i-\-2j =3 

00 

y = -2x 3 + 2/j.xy + bijX l y 3 =Y(x,y). 

i-\-2j=4 


( 2 . 1 ) 


Theorem 2.1. For any positive integer s and a given number sequence {cog}, /3 > 
3, a formal series can be constructed successively in terms of the coefficients c a g 
(a 7 ^ 0) as 

OO OO 

M(x,y) = y 2 + c a px a yP =Y M k(x,y), ( 2 - 2 ) 

a+/3=3 fc=2 

satisfying 


dX 

dx 


dY 

dy 


M- (s + 1) 


fdM 


v dM\ ^ . \ t 

x + ~d^ Y ) = 

* ' m= 3 


(2.3) 


where Mj~{x,y) is a kth-degree homogeneous polynomial in x and y, satisfying 
sp = 0 for all k. 
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Theorem 2.2. For a>l,a+P >3 in (12.21) and (12.31) . c a p can be uniquely 
determined by the recursive formula, 

Cap = 7——TT— (^4ct-l,/3+l + -Ba-l„S+l)- (2-4) 

-L JCK 

For to > 1, w m (s,/x) can fee uniquely determined by the recursive formulae: 


where 


M) — ^4m,0 4“ O 5 

N _ ljJ2m+4(s, fl) 

m ~ 2to - 4s - 1 ' 


a+/3 —1 

A a p = [k - {s + l)(a - k + l)]akjC a -k+i,p-j, 

k+j=2 
a+p — 1 

BaP — ^ [j (^ 4“ 4)(/^ J 4“ l)]fefcjCa — fc,/3 — 4 + l • 

fc+j=2 


(2.5) 

( 2 . 6 ) 


(2.7) 


Theorem 2.3. The origin of system EH is an analytic center if and only if 
the origin of system eh is a center of oo -class, namely, the origin of system 
EH is a center for any natural number s. 


3. Analytic centers of system (11.41) 

Now, we discuss the analytic centers of system EH in two cases. 

3.1. Case 1: n = 1. 

For this case, system EH can be written as 

x = y + a 30 x 3 + a 23 x 2 y + a 32 xy 2 + a 03 y 3 + a 50 x 5 + a 4 ix 4 j/ 

+a 32 x 3 y 2 + a 23 x 2 y 3 + a u xy 4 + a 05 2 / 5 , (3.1) 

y = —2a; 3 . 

According to Theorem 12.11 we can find a formal series M(x,y) = x 4 + y 2 + 
o(( x 2 + y 2 ) 2 ) for system (13.11) . such that (12.31) holds. Applying the recursive 
formulae in Theorem 12.21 to system EH, with the help of Mathematica, we 
obtain 


w 3 — w 4 — W5 — 0, 
w 6 = (4s - l)a 30 , 
w 7 = 3(a + l)c 0 3, 

ui 3 = — — (4s — 3)(2a 4 2 4- 5aso), 

5 

u g = 0 , 

wio = — y(4s — 5)(2a 3 2 + 3a2iaso), 
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Wn = 

Wl2 = 
Wl3 = 
Wl4 = 

Wl5 = 
Wl 6 = 
Wl7 = 
Wis = 

Wig = 

W 20 = 

W 21 = 
W 22 = 

W23 = 

W 24 = 

where 


— (s + l)c 05 , 

-Tt(4s - 7)(12ai4 + 30a O 3a 5 o + 5a 4 ia 50 ), 

45 

0, 

— —^-(4s - 9)(6a 23 + a 2 ia 4 i - 10a 2 0 ), 

35 ,, 

—(S + 1)C 0 7, 

— YYj(4s — ll)(60ao5 + 10ao3a4i + — 3a2ia 2 g), 

0 , 


a 50 

1155 

315 

~64 


(4s - 13)(2a 2 ia4 1 + 300ao3a 2 o + 9a 2 ia 2 0 + 100a4ia 2 o ), 
(s + 1)CQ9, 


«50 

’895050 


(4s — 15)(28a2ia| 1 + 252024044050 + SOOa^agg + 567 ^^^ 

+ 3600a2ia4ia| o + 4500o5 0 ), 


0 , 


^ a5 ° (4s — 17)(4a2 1 a| 1 + 36024044050 + 106021044050 + 81024a 4 


235125 
693 


50 
2 „4 


50 


128 


(s + l)coil, 


Q50 


42089726250000(s+ 1) 


+ 450o2iO4iagg - 125044050 ), 


/1, 


(3.2) 


fi=- 151748682120^041 - 84454927200a2ia|i + 22768748000a^ati - 136573813908a^a2iO§ 0 

- 1193662008000a^atiO^ 0 - 2087643726000a^a|i a 2 50 + 651216400000a 2 ia;(iol? 0 

- 307291081293a^O5 0 - 3661266760200a^a 4 ia^o - 9313504335000a^O4i4) 

+ 5946721200000a^a^a|o + 23826000000a^af 0 - 17785962180a^a|iS 

- 98929404000o^a^s + 26934900000a^a^s - 160073659620a W 50 s 

- 1398534984OOOo2io|iO50S - 2442981870000a^o^a^s + 18810000000o^a^ 0 s 
+ 770106000000a2iatia^s - 360165734145a^O50 s + 6998670000000a!i£&a£ o s 

- 4290086997000a^a 4 ia5oS - 10903637205000a2ia4ia5oS + 4967473392a^a4iS 2 
+ 27628675200a^a^s 2 - 7529168000a^a^s 2 + 44707260528a^ia^ 0 s 2 

+ 390585888000af4atiO^ 0 s 2 + 682203816000af 1 a;^a| ) s 2 - 215262400000o2ia^4a^ 0 s 2 
+ 100591336188a^ 0 s 2 + 1198155823200a^a4ia^ 0 s 2 + 30449730600000^0^0^ s 2 

- 1955419200000a|iatiaf 0 s 2 - 5016000000a|ia| o s 2 . 
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Based on (12.61) and (13.21) . it is easy to find the first ten quasi-Lyapunov constants 
of system m- 

Theorem 3.1. The first ten quasi-Lyapunov constants at the origin of system 
63 are given by 


Ai — 030, 

A2 = -( 2 ai 2 + 5 a 50 ), 

5 

A 3 = -(2032 + 3021050 ), 

A 4 = — (12ai4 + 30003050 + 5044050 ), 

45 

A5 = ' (6023 + 021041 — IOo^q), 

A 6 = — yyy(60oos + 10ao304i + 044 — 3021050 ), 

A7 = — 1155(2021041 + SOOaosoio + 9024050 + 100 a 4 ia^ 0 ), 


As = ~ 895050 ( 28 « 2 i Q 4 i + 252a 21 a 41 a 50 + 800 a|i 0 5 o + 567a^a 50 
+ 3600021041050 + 4500a5o), 

Ao = (4Q21Q41 + 36021044050 + 100021044050 + 81024050 


235125 

+ 450024041050 125044050 ), 


A 10 — — 


a 50 


42089726250000(s+ l)(4s - 19) 


fu 


(3.3) 


where Xk-i = 0 for k = 2, ■ ■ • ,10 have been used in the computation of A 

It follows from Theorem 13.II that the following assertion holds. 

Proposition 3.1. For n = 1, the origin of system JO is an analytic center 
if and only if the following conditions are satisfied: 


030 — O 12 — 032 — 014 — 050 — 0. (3-4) 

Proof. By setting Ai = A 2 = • ■ • = A 10 = 0, it is easy to get the conditions in 
(IQ) . Assume 050 ^ 0, and denote 

/2 =28021044 + 2520240^4050 + 800044050 + 567024050 + 3600021041050 + 4500 o 5 0 , 

fs =4034044 + 36024044050 + 100 a 2 ia 41 a^ o + 81034050 + 450034041050 - 125 o 4 iof 0 - 

(3.5) 


7 







Then, we have 


Ri = Resultant[/ 2 ,/3,a 2 i] 

= 252226880859375a§g(37a| 1 + 36000(44} + 864000(4), 
i ?2 = Resultant[/ 2 ,/i,a 2 i] 

= - 750785873641864353168750000000000000000(4(—879390304066912<4 

+ 47983547106994035360(44) + 49445533255803715842660(4(4 
+ 1456057532744172471928500a| 1 dgo + 6498810664995012399669375a^ 

- 2759277767198304a^s + 169297706825726316960 44o s 

+ 173470743593716632941700o| 1 ag 0 s + 5108374765584631369552500 44) s 
+ 22851124455468570581840625a^s - 20806016094293764s 2 
+ 151804543373289707520 44o s2 + 154422423262500291638820a^a| o s 2 
+ 4547533910484627929569500a^a^s 2 + 204009170845121696548856254s 2 
+ 610343850576768ai?s 3 - 33200851039501326720 44) s3 

- 342149472446615260119004ialo s3 - 1007536381850376947992500a| 1 a^s 3 

- 44967293067883449121031254s 3 + 5738786720578564s 4 

- 49876076993258065920a^ 0 s 4 - 50424620011202182011120a^a^ o s 4 

- 14849734841582645575620004°5os 4 - 6678213700042142946607500a^s 4 

- 2172646904355844s 5 + 18276082133674805760a^a| 0 s 5 

+ 18496906029642804955200 444 + 544720023506226322440000 44s 5 

+ 2448661869754899507450000a^s s + 199146343813124s 6 

- 1701986627481384960 44s 6 - 1721646399680525295360 44s 6 

- 50701299015707667936000 44 s 6 + 2279637270657990507600004s 6 )- 

With the aid of Mathematica, we obtain for Vs € Z + 

G\ = Resultant [i?i, i? 2 , « 4 i] 

= - 182848672642886912449902102931881129741668701171875a^(l + s) 6 (-19 + 4s) 6 
x (12242160594943288477497249258950767957 
+57187190996418911124243473597501985540s 
+84210057837841105190444817587559944702s 2 
+ 22053341878592957414426973876225026580s 3 
-34746447450361087057581921863631440523s 4 

- 7190180552428847800895138514692327280s 5 
+ 8952012886140489676856041653019558112s 6 

- 1982180847477328550724618213150339840s 7 

+138354459536790840295491820367594752s 8 ) 3 ^ 0. 

So there are no solutions for the set of equations, /i = / 2 = /3 = 0, implying 
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that there do not exist other analytic center conditions for system eh if ^50 ^ 

0 . 

Under the conditions in (13.dl) . system (El) becomes 

x = y + a 21 x 2 y + a 03 y 3 + a A1 x 4 y + a 23 x 2 y 3 + a 05 y 5 , 

o (3.6) 

V = —2x 3 . 

Obviously, system (13.61) is symmetric with the y- axis. According to Theorem 11 
in Q, the origin is an analytic center of system (EH). □ 

Proposition 13.11 implies that 

Theorem 3.2. The necessary and sufficient conditions for the origin of system 
EH being an analytic center are determined from vanishing of the first ten 
quasi-Lyapunov constants, that is, the conditions given in Proposition \3.1\ are 
satisfied. 

When the cubic-order nilpotent singular point, 0(0, 0) is a lOth-order weak 
focus, it is easy to show that the perturbed system of m, given by 

x = 5x + y + a 30 x 3 + a 2 \x 2 y + a 32 xy 2 + a 03 y 3 + a 5 ox 5 + a^ix^y 

+ a 32 x 3 y 2 + a 23 x 2 y 3 + a^xy 4 + a 05 J/ 5 , (3.7) 

y = Sy-2x 3 , 


can generate ten limit cycles enclosing an elementary node at the origin of 
system (13791) . 

Theorem 2.2 in [32J implies the following result, 


Theorem 3.3. If the origin of system EH is a 10 th-order weak focus, then 
within a small neighborhood of the origin, for 0 < 5 -C 1, perturbing the co¬ 
efficients of system EH can yield ten small-amplitude limit cycles bifurcating 
from the elementary node 0(0,0). 


Proof. The origin of system (13.71) is a lOth-order weak focus if and only if 


5050 3021050 

030 = 0 ; 012 —-—; 032 =-^—; 

«14 = — — (6003050 + 041050); 

023 = g( —021041 + lOajg); 

005 = — (—10003041 — ah + 3021050); 

a ° 3 = T^TT^T (-20210^ — 9023050 — 100 a 4 iaio). 
toUUO 5 0 J 
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and 


Jn = 


<9(Ai, A 2 , A 3 , A 4 , A 5 , A 6 , A 7 , As, Ag) 

d(a 30, 012, 032 , fll4, «23 j a05, O 03 , 021, O 41 ) 


2048«f 0 


56772622027996875 

„4 3 „4 


(22402304! + 3024a2 1 a| 1 Og 0 + llOOOa^ia^agQ 


+ 13608a2 1 O4 1 O5 0 + 108 9 00a2 1 a4 1 ag O 
+ 20412a^a4ia| 0 + 3118500140^0^ 
+ 200475023010 + 1022625o2ia 4 iaf 0 4 


b 123500a44a| o 
b 783000a2iatia|o 

4500000440^). 


Furthermore, 


Jo 


R 5 = Resultant[/ 2 , — 3 —, « 2 i] 


^50 

= 973840®! - 98391600a^a^ o 


24582976875a|4o|o + 80858250000a^, 

Resultant[i? 5 , i? 4 , 021 ] 

= 1928337060674939567063811915624524516160347438902935330714052761390000 


00000000000000000000000a^o + 0. 


So Theorem 2.2 in [32[ yields the conclusion holds. 


□ 


3.2. Case 2: n> 2. 

For this case, system m can be written as 

x = y + x(a 30 x 3 + a 2 ix 2 y + ai 2 xy 2 + a 03 y 3 + a 2n+3 fiX 2n+3 

+ cL 2n+2 ,iX 2n+2 y + a 2n+ i. 2 x 2n+1 y 2 4-b ai, 2n+2 xy 2n+2 

+ ao, 2 n+ 3 y 2n+3 ) = Xi(x,y), 

y = - 2x 3 . 

Theorem 3.4. Forn > 2, the origin of system (13.81) is at most a (n+4)th-order 
weak focus. If the origin of system (13.81) is a (n + 4)th-order weak focus, then 
within a small neighborhood of the origin, perturbing the coefficients of system 
(13.81) can yield n 4-4 small-amplitude limit cycles around the elementary node 
0 ( 0 , 0 ). 

Proof. For a nilpotent system, in order to study the dynamical behavior in 
the neighborhood of the origin, we could consider u and x 2 to be infinitesimal 
equivalence in the neighborhood of the origin, see [321 ]. Construct a comparison 
system, 

x = y + x(a 2 ix y + a 03 y + a 2n + 2 ,ix y 4-b ai, 2n+2 xy T ) 

= X 2 (x,y), (3.9) 

y = -2a: 3 , 

which shows that the system is symmetric with the x-axis, and so the origin 
0(0, 0) is a center. 
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Next, we compute the determinant of system (1X71) to obtain 


Ji = det 


Xi (x,y) 
X 2 {x,y) 


—2x 3 ' 
—2x 3 


= —2x 4 (a 30 x 2 + a 12 y 2 + a 2n+3t0 x 2n+2 

+a 2n +i, 2 x 2n y 2 H-h a 3i2n x 2 y 2n + ai t2n+2 y 2n+2 ). 


By treating the y and x 2 as infinitesimal equivalence in the neighborhood of the 
origin, we have 


Ji = -2x 4 (a 30 x 2 + a i2 x 4 + a 2n+3 fiX 2n+2 + a 2 „+i, 2^ 2 " +4 
+ • • • + a 3 , 2 n x ‘ in+2 + Ql,2n+2^ 4n+4 ), 


which implies that a 30 , ai 2 , a 2n+3t0 , a 2n +i, 2 , • • • , a 3j2n , ai, 2 n +2 could be taken 
as the focus values of system m- So for n > 2, the origin of system (13.81) 
is at most an (n+4)th-order weak focus. According to Theorem 4.1.5 in [311 ] . 
within a small neighborhood of the origin, perturbing the coefficients of system 
(13.81) can yield n + 4 small-amplitude limit cycles around the elementary node 
0 ( 0 , 0 ). □ 
Furthermore, similar to Proposition 13.II we have the following result. 

Theorem 3.5. For n > 2, the origin of system (13.81) is an analytic center if 
and only if 


030 — Ol2 — O2n+3,0 — Cl2ra+1,2 — • • • — d 3i 2n — dl,2n+2 — 0. 


(3.11) 


Proof. When a 3 o — 012 — d 2 n+ 3 ,o — « 2 ti+i ,2 — ■ • ■ — a 3j2n — ai t2n + 2 — 0, 
system m could be rewritten as 

x = y + a 21 x 2 y + a 03 y 3 + a A1 x A y + a 23 x 2 y 3 + a 05 y 5 
+ • • • + a 2 n+ 2 ,lX 2n+2 y + • • • + d0,2n+32/ 2n+3 ) (3-12) 

y = —2x 3 . 

Obviously, system (13.121) is symmetric with the y- axis. According to Theorem 
11 in [)||, the origin is an analytic center of system (13.81) . □ 


4. Analytic centers of system (11.51) 

Now we turn to discuss the analytic center conditions for system (11.51) . It also 
has two cases. 

4-1. Case A: n = 1. 

For this case, system m can be written as 


x = y, 

y = -2x 3 + a 2 \x 2 y + a\ 2 xy 2 + a 03 y 3 + a 50 x 5 + a 4 ix 4 ?/ (4.1) 

+ a 32 x 3 y 2 + a 23 x 2 y 3 + a 14: xy 4 + a 05 y 5 , 


11 






for which we can find a formal series M(x, y) = x 4 + y 2 + o((x 2 + y 2 ) 2 ) according 
to Theorem 12. 11 provided that (12.31) holds. Carrying out calculations with help 
of Mathematica and applying the recursive formulae in Theorem 12.21 to system 

031), we obtain 

W3 = W4 = W5 = 0, 

We = "(4s - 1 )021, 

UJ 7 = 3 (s + l)c 03 , 

w 8 = — i(4s - 3)(6a 03 + a 4 i), 

UJ 9 = 0 , 

wio = — y(4s — 5)(2ao3ai2 — 2023 + 3003050 ), 

15 , 

Wli — — (s + l) c 05, 

x>i 2 = — (4s — 5)(40ao5 — 4003032 — 20 , 030120.50 — 5003050 ), 

W 13 = 0 , 

W 14 = —y(4s — 9)(48ao3 — 40ai4 + 12012032 + 6Oi 2 O50 
154 

+ 12032050 + 21012050 + 18a| 0 ). 

Then, for 012 + 2aso ^ 0, 

35 ^ 

<*>15 = -g-(s + l)c 0 7 , 

<*>16 = ^ 7 ^ (4s — ll)(64oQ 3 ai2 + 16032 + 128oq 3 O50 + 16012032050 


W 17 = 0 , 

wi8 = — 


4012O50 + 32032050 + 20012050 + 23a| 0 ), 

(4s — 13)(4a 3 2 + 2012050 + 5a 2 0 )(ll2a 2 2 ci32 ~ 432a 32 


61600(ai2 + 2050) 

+ 56a 32 O50 — 96012032050 + 200o 2 2 a^ 0 — 640 o32O 2 q 
+ 120ai2Og 0 - 85ag 0 ), 


315 

<*>19 — ~ FT U + s ) c 09, 

64 


W 20 — — 


ao3 


r (4s — 15)(4a 32 + 2012050 + 5al 0 )(14:372996a 4 2 a 3 2 


40840800(ai2 + 2a 50 ) 2 

- 63894256a? 2 032 + 34076160 a \ 2 + 7186498o^ 2 a 5 o - 107341160^032050 

- 12772032oi203 2 a 5 o + 28572751at 2 a5 0 - 99036264ai 2 a32a5o 

+ 45544768a 32 O50 + 26 9 58 1 96a? 2 ai 0 - 39087216oi2a 3 2af 0 ), 


W 21 = 0, 


(4.2) 
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W22 - 11639628000°( 3 ai 2 + 2a 5O ) 2(4a32 + 2012050 + 5a ' o) U ' 
and for ai 2 + 2aso = 0, 

wi6 = |^(4s - 11)(-4a 3 2 + aJo)(4a 3 2 + a§ 0 ), 

W 17 = 0 , 


and in addition if <232 = 

9 ao 3 , 


4 > 


wi 8 — ^^(4s ~ 13)(4a 32 + a 2 o)(24aQ 3 + af 0 ), 
315 . 

W 19 = — (s + !)C09, 

7 a 03 /. 1 r\ 6 

w 20 = -TTT7T7^(4 s - 15)o 50 , 


13260' 


W 21 = 0, 
0 J 22 = — 


003^50 


13856700(1 + s) 


(5391 - 205861s + 66718s 2 ); 


if a 32 = 


where 


315, 

Wl9 — 77 l s + l)Cog, 
64 


W18 = 0 , 

31 
"6 
_ 2 ao 3 
W20 ~ 5525 
w 2 i = 0, 

W 22 = 


(4s - 15)(16aQ3 + a^ 0 )(27a§ 3 + 2a| 0 ), 


4 003050 


1154725(1+ s) 


(16ag 3 + 03 q)(270q3 + 2050 ), 


fi = - 175151884140096a? 2 a 3 2 + 900479057104608a? 2 a^ 2 - 870691837997952oi 2 a| 2 

- 875759420700480^050 + 191710117401504ai 2 a32a 5 o + 108696440458488a 2 2 03 2 a 5 o 

- 135343772601984a| 2 a 5 o - 348204560750520af 2 a^ o + 1449441463187484a? 2 a 32 al? 0 

- 1227387674864544ai2a^ 2 a5 0 - 328545834076764a? 2 af 0 + 705733010157654a? 2 a 32 a! o 

- 180965555611680a 2 2 af o - 196530256579516a? 2 a 32 s + 956545497558256a? 2 a^ 2 s 

- 775708121404800ai 2 a| 2 s - 98265128289758a? 2 a 5 os + 188198718433828at 2 032a 5 os 

- 92031891970176a| 2 a 5 os - 390699835897045af 2 a 2 0 s + 1519172249991516a? 2 a32a 2 0 s 

- 1080068393470816a 1 2a 2 2 o 2 0 s - 368634374118786af 2 oi 0 s + 690504693129726a? 2 a 3 2aios 

- 123060065895200o3 2 Og 0 s + 80500800862640a^ 2 a 3 2S 2 - 396831061224512af 2 a^ 2 s 2 

+ 336523546570752oi 2 o| 2 s 2 + 40250400431320a? 2 a 5 oS 2 - 79597335692936aV32a 50 s 2 

- 56511998888512a 2 2 a 2 2 a 5 os 2 + 43311880631808a| 2 a 5 os 2 + 160035098537960af 2 a^ o s 2 

- 632261525950008af 2 a 32 a 2 o s 2 + 470151299998208ai 2 a| 2 o| 0 s 2 + 150997488382038at 2 ai o s 2 
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+ 142715605813496ai 2 a3 2 a 5 os - 292287689531688oi 2 a 32 ai?oS 2 + 57905489716480a 2 2 af 0 s 2 . 

Based on (USD and AMD , it is easy to find the first nine quasi-Lyapunov constants 
of system (ED- 

Theorem 4.1. The first nine quasi-Lyapunov constants evaluated at origin of 
system G 3 D are given by 

Ai = —-o 2 i, 

A 2 = —-(6ao3 + 041), 

5 

A 3 = —-(2ao3ai 2 — 2a 2 3 + 3003050 ), 

A 4 = — (40ao5 — 4ao303 2 — 2ao3ai 2 a5o — 5 003050 ), 

A5 = (48oqs — 40ai4 + 12ai 2 a 32 + 6a 2 2 a5o + 12a3 2 a5o + 21ai 2 a 2 0 + 18053). 

Then, for ai 2 + 2 aso P 0 , 

CL 03 / O q n q q q 

A 6 = ■^j(64a 03 ai2 + 16a 32 4- 128a 03 a5o + 1^(112^32^30 + 4a 12 a 50 + 32a32a 50 

+ 20 a 12 a| 0 + 23a5 0 ), 

A 7 = — g „ nn , ———~-r(4a3 2 + 2ai 2 a5o + 5a 2 o)(112a 2 2 a3 2 — 432a 2 2 + 56ai 2 a5o 

61600(ai2 + 2050) 

— 96012032050 + 200a 12 a 50 — 640032053 + 120012053 — 85053), 

As = ~ 40840800(ai 2 + 2os 0 ) 2 (4 ° 32 + 2aM + 5<4)(14372996of 2 « 32 

- 63894256ai 2 a 2 2 + 34076160a| 2 + 7186498af 2 a 5 o - 10734116a? 2 a 3 2a5o 

— 12772032ai 2 a 2 2 a5o -I- 28572751 oi 2 o 2 o — 99O36264o 2 2 O3 2 o 2 0 
+ 45544768a 2 2 a5o + 26958196a? 2 a| 0 - 39087216ai 2 a 32 afo), 

A9 = 11639628000°(a 12 + 2o 50 ) 2 (4 ° 32 + 2ai2aso + 5a '°) ^ 


while for 012 + 2aso = 0, 


As — 


Q 03 

520 


(—4a 32 + a 2 0 )(4a32 + a 2 0 ), 


and in addition if 032 = pp, 


A 7 

As 

Ag 


9ao3 

7700' 


(4a 32 + a 2 o)(24ao3 + 053), 


7 Qq3 
13260 


a 


6 

50 1 


O03Q50 

13856700(1 + s) 


(5391 


205861s + 66718s 2 ); 
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*/a 3 2 = — f-, 


Ar = 0, 

2 ao3 


As — 


5525 


An — —- 


(16a 03 + a§ 0 )(27a 03 + 2a\ 0 ), 
4 003050 


1154725(1 + S ) (16 °- + a so )(27a 03 + 2 ( 4); 

where \k-i = 0 for k = 2, • • ■ ,9 /icroe been used in computing A 
Furthermore, the following result can be easily obtained. 


Proposition 4.1. For n = 1, origin of system G 3 D is an analytic center if and 
only if one of the following conditions holds: 



«21 — 0-03 — a 41 — 023 

— 0 05 — 0; 

( 4 - 3 ) 

a 2 i 

= 014 = ao 5 = 0, a 4 i 

= —6003, 

( 4 - 4 ) 

a 2 3 

= 5(2012 + 3050)003, 

„2 __ 1_ 3 . 

U 03 — I 6 a 50 i 

021 =ao5 = 

0, a 41 = —6003, o 12 = - 

-2050 , 0 2 3 = —5003050, 


1 2 

2 2 3 

1 3 

( 4 . 5 ) 

a 32 — 4 a 50 > a 03 ~ 27° 50 ’ ai4 — 

72 « 5 °- 



Proof. It is easy to get the conditions (11.3D - (14.5D by setting Ai = A 2 = • • • = 
A 9 = 0. When 050 0, let 

/s = 112a? 2 a 32 — 432a 32 + 56a^ 2 fl50 — 96(112032050 + 200a 2 2 a5Q — 640a 32 O5Q 

+ 120ai 2 a5Q — 85a5g, 

/ 6 = 14372996af 2 a 32 - 63894256a? 2 a§ 2 + 34076160(4 + 7186498a^ 2 a 50 

- 10734116a 32 a 32 a5o - 12772032ai 2 a^a 5 o + 28572751(44 

- 99036264a? 2 a 32 a^ o + 45544768a§ 2 a§ 0 
+ 26958196a 32 a5Q — 39087216ai2a 32 a5 O . 

Then, we obtain 


R 3 = Resultant [/ 4 , /s, a 32 ] 

= - 200206652313600a^ o (ai2 + 2a 50 ) i {72a 5 12 + 652a? 2 a 50 + 2694a? 2 a5 0 
+ 6043 a 2 12 a 3 50 + 7092ai 2 af o + 3463a| 0 ), 
i? 4 = Resultant[/ 4 , f$, a 32 ] 

= - 289298612593152000af o (ai2 + 2a 5O ) 4 (278734084992aj' 2 + 3588989330016a? 2 a 50 
+ 21026493958464a? 2 4 + 71854303647672af 2 a| 0 + 152324731255716a? 2 a? 0 
+ 197669760539760a? 2 a^ o + 144211704399495ai 2 af o + 45502270176438a£ o 
+ 1651136278704a 32 s + 18838630495800a? 2 a 5O s + 98153819816532a? 2 a§ 0 s 
+ 294846821571666a? 2 aj? 0 s + 534120897148782a? 2 a? o s + 565346904516561a? 2 al( o s 
- 310952885702031ai 2 af o s + 62534145621954a£ 0 s + 1859553268056aj’ 2 s 2 
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+ 17577034105268a? 2 a 50 s 2 + 76054378966082a'^ o s 2 + 181916843290105af 2 af 0 s 2 
+ 238094684972342a? 2 a 4 o s 2 + 147267056136244a 2 2 at 0 s 2 + 19507558179230ai 2 a^ o s 2 

- 6987410240074aJ 0 s 2 - 1714957345728a]> 3 - 17270801713568a? 2 a 5O s 3 

- 80278393721648a 5 2 a 2 0 s 3 - 212296048411048a([ 2 af 0 s 3 - 328607615124608a? 2 a^ o s 3 

- 285216632755972a 2 2 a| 0 s 3 - 121914622605938ai 2 af 0 s 3 - 19568176640258a£ o s 3 
+ 301352777088a[ 2 s 4 + 3072037299008af 2 a 50 s 4 + 14465726297408af 2 a^ 0 s 4 

+ 38896515312256af 2 af 0 s 4 + 61684793716376a? 2 a^ 0 s 4 + 55602984609100a? 2 af 0 s 4 
+ 25319000517368ai 2 af 0 s 4 + 4451109045332a^ o s 4 ). 

Further, with the aid of Mathcmatica, we obtain for Vs £ Z + 

G 2 = Resultant [-R 3 , i? 4 , ai 2 ] 

= 30984189289342953910272000a 3 o(l + s) 5 (-17 + 4s ) 5 
x (-123287750793562256929839075859953216 
-33902812452688795016920021129342044624s 

- 180855325034978657368019444342423080236s 2 

- 1067066959204615961659004741488392865575s 3 

- 3328343437962444375340762099992891472110s 4 

- 4773196954655562390848005555854946921459s s 
+14241540803784759916727236436410335714320s 6 

- 10732088467496096467063502795815305475120s 7 
+ 3721436248399857364295558363131840668032s 8 
-625676462230475741935920059840273863168s 9 

+ 41454785818979861302571809000901003264s 10 ) ^ 0. 

The above calculations indicate that the equations /4 = /s = /6 = 0 do not 
have real solutions, namely, there do not exist other analytic center conditions 
for system (ED If ^50 7 ^ 0 . 

When the conditions in (14.31) hold, system (14.11) becomes 

± = V > (4.6) 

y = -2x 3 + ai 2 xy 2 + a 50 x 5 + aj, 2 x 3 y 2 + a^xy 4 . 

Obviously, this system is symmetric with the y-axis, implying that the origin of 
system (14.61) is an analytic center due to Theorem 11 in [9j. 

When the conditions in (14.41) are satisfied, system (14.11) becomes 

x = y, 

y = ^(- 8 :r 3 + 4a 50 a ; 5 - 24a 03 a; 4 y + 4a i 2 a;y 2 - 2a 12 a 50 x 3 y 2 (4.7) 

- 5a 2 0 x 3 y 2 + 4a 03 y 3 + 4 a 03 ai 2 a; 2 2/ 3 + 6 a 03 a 50 a; 2 y 3 ). 
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Introducing the transformation, 


V = 


and time scaling, 


T = 


(-2 + a 50 x 2 )z 
2(—2 + a 5 ox 2 + ao 3 Xz) ’ 

2(—2 + a 5 oX 2 ) 3 t 
-2 + CL 50 X 2 - 2ao3 xy ’ 


into system (Id.71) results in 
dx 


— = z{a 5 0 x 2 - 2) 2 , 
dz 1 

— = — -x(l28x 2 — 192asoX 4 + 96a 2 0 x 6 — 16af 0 x 8 — 16ai 2 2; 2 


(4.8) 


+ 24ai2a5ox 2 z 2 + 20 a 2 o x 2 2 : 2 — 96ao 3 x 4 ^ 2 — 12 ai 2 a 2 0 x 4 2 : 2 
— 20a 3 0 x 4 £ 2 + 48ag 3 a5ox 6 2; 2 + 2ai2a 3 0 x 6 2: 2 + 5a 4 0 x 6 z 2 
+ 32ao 3 x 5 z 3 + 2 a 0 3 a| 0 x 5 z 3 ), 

which is symmetric with the z-axis because Oq 3 = — Thus, according to 
Theorem 11 in [ 9 }, the origin of system (1X71) is an analytic center. 

Similarly, when the conditions in (14.51) hold, system (14.11) becomes 


32 


x = y, 

y = 144x 3 + 72a 50 x 5 - 432a 03 x 4 y - 144a 50 xy 2 - 18a 2 0 x 3 y : 

+ 72 a 03 y 3 - 36a 03 a 50 x 2 y 3 - a 3 0 xy 4 ), 
for which there exists an analytic integrating factor, 


(4.9) 


i(x,y) = 


- In 2 X 4 
£ 8 U 50 J ' 


(1 - \a 50 x 2 + |a 50 xy) 4 ’ 


indicating that the origin of system m is an analytic center. 
Therefore, Proposition 14.11 implies the following result. 


□ 


Theorem 4.2. The necessary and sufficient conditions for the origin of system 
HD being an analytic center are given by the vanishing of the first nine quasi- 
Lyapunov constants, that is, one of the conditions in Proposition^. l\is satisfied. 


Similarly, when the cubic-order nilpotent singular point 0(0, 0) is a 9th-order 
weak focus, it is easy to prove that the perturbed system of hid, given by 


x = Sx + y, 

y = Sy - 2x 3 + a 2 ix 2 y + ai 2 xy 2 + a 03 y 3 + a 50 x 5 + aux A y 
+ a 32 x 3 y 2 + a 23 x 2 y 3 + ai 4 xy 4 + a 05 y 5 , 


(4.10) 


can generate nine limit cycles enclosing an elementary node at the origin. The 
proof is similar to that for Theorem 6. 
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Theorem 4.3. If the origin of system (14.101) is a 9th-order weak focus, then 
within a small neighborhood, of the origin, for 0 < 6 <C 1, system (14.101) can 
yield nine small-amplitude limit cycles around the elementary node 0(0,0). 

Proof. The origin of system (14.101) is a 9th-order weak focus if and only if 

021 = 0 ; 

041 = ~6ao3; 

^23 = 2°03(2ai2 + 3 o5o); 

a 05 = -^(^ a 03 a 32 + 2 003042050 + 5003050); 

3 

Ol4 = — (16003 + 4012032 + 2 Oi 2 O50 + 4032 050 + 70120^ + 6050); 

°03 = Q4(a 4 - 2 a )^~ — 16012032050 — 4a 12 O50 — 32o32a 5 o — 20012050 — 23a 5 o). 

and 

j _ d(Xi, A 2 , A 3 , A 4 , A 5 , A 6 , A 7 , As) 

<9(a 2 l, 041, a 2 3, 005, Oi4, Oo3, 032, Oi 2 ) 

= ~ 38846808000(012 + So) 3 (16 ° 32 + 16ai2a32 ° 50 + + 32a32 °'° + 20ai2a '° + 23a ^ 

x (9658653312af2a 2 2 - 94084036480af 2 a| 2 + 233154747648a? 2 a^ 2 

- 58883604480o3 2 + 9658653312ai 2 O3 2 a50 — 65797494336af 2 o| 2 a5o + 14579296576af 2 O3 2 O50 
+ 288547550208ai 2 a^O50 + 2414663328a? 2 a§ 0 + 28912166304af 2 a 3 2a^o 

- 406316235904af 2 a| 2 aj? o + 936674204416a? 2 a| 2 a| ) - 127467869184a^ 0 
+ 19144952176a( 2 a| 0 - 109760863872af 2 a 3 2afo - 244778217136a? 2 a 2 2 af 0 
+ 1036308069632ai 2 a| 2 aio + 49040513184a? 2 o| 0 - 522094417624af 2 a 3 2af 0 
+ 985367488400a 2 2 a 2 2 O50 - 6731934720al 2 O5 0 + 18296999196af 2 a^ 0 

- 545453749756a? 2 a 3 2a|o + 1124651917440ai 2 a| 2 a|o - 114440868400af 2 a^ o 
+ 131607099120a? 2 a 3 2a^o + 107989619520a 2 2 af o - 173404183515a? 2 a^ o 

+ 352746916240ai 2 a 3 2o£o - 62864473110a? 2 a| o + 39921858880a 32 afo + 12556768140ai 2 al o )- 

Furthermore, 

i ? 6 = Resultant[/ 5 , J 2 ,a 3 2 ] 

= (49a 3 2 + 196oi 2 O50 + 35 7ai2O5 0 + 466og 0 ) 2 x (225792 oi 2 + 3214400ai 2 a5o 
+ 21345968o( 2 a 2 o + 87227792 a 6 12 a 3 50 + 246151320af 2 a5 0 + 509498952ai 2 a5o 
+ 785491407a 3 2 a^o + 864775852a? 2 a^ 0 + 601885864ai 2 a| o + 196001002a^ o ), 

Resultant [Rq, R 3 , 012 ] 

= - 110020282897692638814502001588125819626479649518678944188134258621907 
66734966784000000a5o^ 0. 
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So Theorem 2.2 in [32[ yields the conclusion holds. 


□ 


f.2. Case B: n > 2. 

For this case, system (USD can be written as 
x = y, 

y = -2x 3 + {a 21 x 2 y + a 12 xy 2 + a 03 y 3 + a 2n+3fi x 2n+3 

+ a 2n+2 ^x 2n+2 y + a 2n+1 , 2 x 2n+1 y 2 H-b a li2n+2 xy 2n+2 (4.11) 

+ U0,2n+3y 2n+3 ) 

= Yi(x,y). 

Theorem 4.4. Forn > 2, the origin of system (14.111) is at most a (n-\A)th- order 
weak focus. If the origin of system (14.111) is a (n+ 4 ) th-order weak focus, then 
within a small neighborhood of the origin of its perturbed system,, perturbing 
the coefficients of system (14.111) can yield n + 4 small-amplitude limit cycles 
enclosing the elementary node 0(0,0). 

Proof. The proof is similar to that for Theorem l3.41 We construct a comparison 
system for system OHD, 


x = y, 

y = -2x 3 + x(a 12 xy 2 + a 2n+3y0 x 2n+3 4-f a 2t2n+ ix 2 y 2n+1 ) (4.12) 

= Y 2 {x,y). 

It is easy to see that system (14.121) is symmetric with the x-axis, and so 0(0, 0) 
is a center. 

Next, we compute the determinant of system (14.121) . yielding 


J 3 = det 


y Yi{x,y) 
y Y 2 {x,y) 


= a 21 x 2 y 2 + a 03 y 4 + a 2n+2 ^x 2n+2 y 2 + a 2rit3 x 2n y i 
H-b a 2t2n+ ix 2 y 2n+2 + a 0y2n+3 y 2n+4 . 


Similarly, we take the y and x 2 as infinitesimal equivalence in the neighborhood 
of the origin in order to study the dynamical behavior of (14.111) around the 
origin. So, J 2 becomes 


J 3 =x 4 (a 23 x 2 + a 03 x 4 + a 2n+2 jx 2n+2 + a 2n , 3 x 2n+2 
+ • • • + a 2i2n +ix 4n+2 + a 0 , 2 n+ 3 x 4 ” +4 ), 


(4.13) 


implying that a 2 i, a 03 , a 2n+2 ,i, a 2n>3 , • • • , a 2j2n+ i, a 0)2 „ + 3 could be considered 
as the focal values of the system. Therefore, for n > 2, the origin of system 
(14.111) is at most a (n+4)th-order weak focus. According to Theorem 2.2 in [321 ]. 
within a small neighborhood of the origin, one can perturb the coefficients of 
system (14.111) to obtain n+4 small-amplitude limit cycles around the elementary 
node 0(0,0). □ 

Moreover, we have a similar theorem for this case. 
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Theorem 4.5. For n > 2, the origin of system (14.111) is an analytic center if 
and only if 


021 — 003 — 02n+2,l — <22n,3 — • • • — 02,2n+l — O0,2n+3 — 0. (4-14) 

Proof. When a ,21 = O 03 = 02 n+ 2 ,l = Q2n,3 = ■ • ■ = < 22 , 2 n+l = Oo,2n+3 = 0, 
system (14.111) can be rewritten as 

x = y, 

y = ~2x 3 + a 12 xy 2 + a 2n+3fi x 2n+3 + a 2n+1 , 2 x 2n+1 y 2 H-f a lt2n+2 xy 2n+2 , 

(4-15) 

Obviously, system (14.151) is symmetric with the y- axis. According to Theorem 
11 in Q, the origin is an analytic center of system (14.111) . □ 

5. Conclusion 

In this paper, two classes of lopsided systems have been studied on their analytic 
integrable conditions and bifurcation of limit cycles. We have obtained some 
analytic integrability conditions for each class of the systems for case n = 1. 

By using certain transformations or integrating factors, we have proved that all 
conditions are sufficient and necessary. For case n > 2, we have constructed 
different comparison systems for each class of the systems and shown that n + 4 
limit cycles may bifurcate from the origin of each system. In addition, conditions 
for the origin being an analytic center are obtained simultaneously. 

6. Appendix 

Detailed recursive MATHEMATICA code to compute the quasi-Lyapunov 
constants at the origin of system (13): c [0,0]=0, c [1, 0]=0, c [0, 1]=0, c [2, 

0]=0, c [1, 1]=0, c [0, 2]=1; when kjO or jjO, c [k,j]=0; else 

c[k, j] = y(-10a 5O c[-5 + k,l+j] + a 50 kc [-5 + k,l+j]~ 5a 50 sc[-5 + fc, 1 + j] 

+ a5oksc [—5 + k, 1 + j] — 8a4ic[—4 + fc, j] + ankc [—4 + k, j] — 4a4isc[—4 + fc, j] 

+ a 4 ifcsc[—4 + k, j] — 4c [—4 + k, 2 + j] — 2jc[—4 + k,2 + j] — 4sc[—4 + k, 2 + j] 

— 2jsc[—4 + k, 2 + j] — 6a 32 c[—3 + fc, —1 + j] + a 32 kc [—3 + k, — 1 + j] 

— 3a 32 sc[—3 + k, — 1 + j} + a 32 ksc [—3 + k, — 1 + j} + a 30 ksc [—3 + k, 1 + j] 

— 6a 30 c[—3 + k,l + j] + a 30 kc [—3 + k,l + j] — 3a 30 sc[-3 + k, 1+ j} 

— 4tt23c[—2 -)- k : —2 -f- j] + a 23 kc [—2 + k , —2 + j] 2 u 23 ^c[—2 -j- k : —2 -)- j] 

+ a 23 ksc[—2 + k, — 2 + j] - 4a 2 ic[—2 + k, j] + a 2 \kc[—2 + k,j] - 2a 2 isc[—2 + k, j] 
+ a 2 \ksc[—2 + k,j] — 2ai4c[—1 + k , —3 + j] + ai 4 fcc[—1 + k, — 3 + j] 

— ai 4 sc[—1 + k, — 3 + j} + auksc[— 1 + k, -3 + j ] — 2ai 2 c[—1 + k, — 1 + j] 

+ ai 2 kc [—1 + fc, — 1 + j] — ai 2 sc[-l + k.—l+ j] + a\ 2 ksc[— 1 + k,—l + j] 

+ a 05 kc[k , —4 + j] + a 05 ksc[k , -4 + j] + a 03 kc[k , -2 + j] + a 03 ksc[k.. -2 + j]). 
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w[m\ = 9a5oc[—4 + to, 0] — asoTOc[—4 + to, 0] + 4a5oSc[—4 + to, 0] — asoTOsc[—4 + to, 0] 

+ 7 a 4 ic[—3 + m, —1] — <Z 4 iTOc[—3 + m, —1] + 3 o 4 isc[—3 + m, —1] 

— a 4 imsc[—3 + m, — 1 ] + 2c[—3 + to, 1 ] + 2sc[—3 + to, 1 ] + 5a32c[—2 + m, —2] 

— a 32 mc[—2 + m, —2] + 2 a 32 sc[—2 + to, —2] — a 32 msc[-2 + m, —2] 

+ 5a3oc [—2 + to, 0 ] — a 3 omc [—2 + to, 0 ] + 2 a 3 osc [—2 + m, 0 ] — a 3 oTOSc [—2 + m, 0 ] 

+ 3a23c[—1 + to, —3] — a 23 TOc[-l + m, —3] + a. 23 sc[— 1 + m, —3] — a 23 TOSc[-l + to, —3] 

+ 3a2ic[—1 + TO, —1] — a 2 lTOc[—1 + TO, —1] + <l 2 lSc[—1 + TO, —1] — a 2 lTOSc[—1 + TO, —1] 

+ oi 4 c[to, —4] — ai 4 rnc[m, —4] — oi 4 TOSc[to, —4] + ai 2 c[ro, —2] — ai 2 toc[to, —2] 

— ai 2 rnsc[TO, —2] — ao 5 c[l + to, —5] — ao 5 mc[l + to, —5] — aossc[l + to, —5] 

— ao 5 TOSc[l + m, —5] — ao 3 c[l + to, —3] — ao 3 TOc[l + m, —3] — ao 3 Sc[l + m, —3] 

— ao 3 TOSc[l + TO, —3] — C[1 + TO, —1] — TOc[l + TO, —1] — Sc[l + TO, —1] — TOSc[l + TO, —1]. 

Detailed recursive MATHEMATICA code to compute the quasi-Lyapunov 
constants at the origin of system (24): c [0,0]=0, c [1, 0]=0, c [0, 1]=0, c [2, 

0]=0, c [1, 1]=0, c [0, 2]=1; when kjO or jjO, c [k,j]=0; else 

c[k,j} = - ^ ^ ( 2 a 50 c [-6 + fc, 2 + j] + a 50 jc[-Q + k, 2 + j] + 2 a 50 sc [-6 + fc, 2 + j] 

+ a,5ojsc [—6 + fc, 2 + j] + a 4 ijc[—5 + fc, 1 + j] + a 4 isc[—5 + fc, 1 + j] 

+ « 4 ij'sc [-5 + fc, 1 + j] - 2a 32 c[—4 + k,j) + a 32 jc [-4 + k,j] + a 32 jsc [-4 + k,j ] 

— 4c[—4 + fc, 2 + j] — 2jc [—4 + fc, 2 + j] — 4sc[—4 + fc, 2 + j] — 2jsc [—4 + fc, 2 + j} 

— 4 a 2 3 c[—3 + fc, -1 + j] + a 23 jc[-3 + fc, -1 + j] - a 23 sc[-3 + fc, -1 + j] 

+ a 23 jsc [-3 + fc, -1 + j] + a 2 ijc[—3 + fc, 1 + j) + a 2 isc[-3 + fc, 1 + j] 

+ « 2 i jsc[—3 + fc, 1 + j] - 6 ai 4 c [—2 + fc, -2 + j} + aujc[—2 + fc, — 2+ j} 

— 2 ai 4 sc [—2 + k,—2 + j] + aujsc[—2 + fc, — 2 + j] — 2 ai 2 c [—2 + fc, j] + a 12 jc [—2 + fc, j] 
+ ai 2 jsc[-2 + fc, j] - 8 a 0 5 c[-1 + fc, -3 + j] + a 05 jc [-1 + fc, -3 + j] 

— 3 ao 5 Sc[-l + fc, -3 + j] + ao 5 jsc[—1 + fc, -3 + j] - 4ao 3 c[-l + fc,-1 + j] 

+ O 03 jc[-l + fc, -1 + j] - a 03 sc[-l + fc, -1 + j] + a 03 jsc[-l + fc, -1 + j]). 

w[m] = — O 5 oc[—5 + TO, 1] — O 50 Sc[—5 + TO, 1] + a 4 ic[—4 + TO, 0] + 3a32c[—3 + TO, —1] 

+ a 32 Sc[—3 + TO, —1] + 2c[—3 + TO, 1] + 2sc[—3 + TO, 1] + 5a23c[—2 + m, —2] 

+ 2a2 3 sc[—2 + to, —2] + a 2 ic[—2 + to, 0] + 7ai4c[—1 + to, —3] + 3ai4sc[—1 + to, —3] 

+ 3 ai 2 c[—1 + m, —1] + a 12 sc[— 1 + m, —1] + 9 a 0 5 c[m, —4] + 4 a 0 5 sc[m, -4] 

+ 5oo 3 c[to, —2] + 2 ao 3 Sc[m, —2] — c[l + to, —1] — mc[l + m, —1] 

— sc[l + m, — 1 ] — msc[l + to, — 1 ]. 
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